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Abstract This article makes use of the apparent indifference that the market 
has been devoting to the developments made on the fundamentals of quantita- 
tive finance, to introduce novel insight for better understanding market evolution. 
We show how these drops and crises emerge as a natural result of local economi- 
cal principles ruling trades between economical agents and present evidence that 
heavy-tails of the return distributions are bounded by constraints associated with 
the topology of agent relations. Finally, we discuss how these constraints may be 
helpful for properly evaluate model risk. 
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1 Introduction 

Most of the research made on financial and economic variables has been dom- 
inated by an epistemological approach, in which the behavior of the economic 
system is explained by a few key characteristics of the behavior itself, like the 
amplitude of price fluctuations or the analytical form of the he avy-tailed return 



amplitude 01 price fluctuations or tne analytical torm 01 tne ne avy-tailed returi 
distributions CMantegn a and StanlevLfl995HBorland et al.Ll2005h . These key char 



acteristics motivated researchers to assume the distribution of the return of those 
variables as a-stable Levy dis tributions or truncated a-stable Levy distributions 
( Manteena and Stanley! . |l994h . The logical reason for this assumption is given by 



the Central Limit Theorem that states that the aggregation of a growing number 
of random variables converges to an a-stable Levy distribution. If these random 
variables have finite variances then the resulting aggregation is a 2-stable Levy 
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distribution, i.e. a Gaussian distribution. If the variance of the random variables 
are infinite, or in the system size order, then the resulting distribution is a 2-stable 
Levy distribution and the so called heavy-tailed distributions appear as a result 
of the aggregation. Further, from an epistemological point of view it is reasonable 
to assume that measurements on aggregates of human activities will result in a 
a— stable Levy distribution, since humans are too correlated with each other to 
show Gaussian "behavior". Henceforth, we refer to a-stable Levy distributions 
with a < 2 as Levy distributions and with a = 2 as Gaussian distributions. 

Without leaving an epistemological approach, we could address the study of 
the resulting distributions by ignoring the previous arguments and construct a 
function that fits any set of empirical data just by building up fitting parame- 
ters until the plotted function fit the empirical data. Such approach would be 
the best one, if economic processes were stat i onary. Unfortunately the y are not 
(|Bouchaud and Pottersl 120031 ISornettel 120031 Ida Cruz and Lindl |2010|) and this 
means that we cannot disregard the underlying mechanisms in the analysis of 
empirical data. Consequently, returning to the above arguments, we are mathe- 
matically limited to a-stable Levy distributions. 

Since heavy-tails are observed in the returns of economic variables, one would 
expect that practitioners use Levy distributions. Why do they use Gaussian distri- 
butions instead? There are good reasons for that. Assuming that Levy distributions 
are good representations of economic variables fluctuations, a model based on it is 
closed when one fits the distribution to empirical data to obtain the remaining con- 
stants, to turn the mathematical model into a financial one. However, a a-stable 
Levy distribution is not necessarily a good representation of future fluctuations. 
There is no guarantee that today's fitting will be the same as tomorrow's. There- 
fore, since working with a Gaussian curve is more straightforward than working 
with a Levy distribution and needs less parameters for curve fitting, there is no 
practical gain in abandoning Gaussian distribution to model the distribution of 
fluctuations according to a prescribed mathematical model, even though it is not 
entirely correct. That is, if a Levy distribution is fitted to empirical data of a non- 
stationary process one will carry basically the same model risk, as if a Gaussian 
distribution is used. 

The differences between Gaussian and Levy distributions as possible solutions 

for the distribution of fina ncial variable values find their fundamentals in the Cen- 

tral Limit Theorem (CLT) (jEmbrechts et al.l.ll997llBillingslevlll995llMantegna and Stanley!. 
1995). In its classical and stronger version CLT states that a sequence of indepen- 
dent and identically distributed random variables with finite expected values and 
variance, has its sample average approaching a Gaussian distribu tion when the 
sample size is increased. Lyapunov version proves (|Billingsle vl ll995h that the same 
is valid in the case where the random variables are not identically distributed, re- 
quiring that they have moments of order greater than tw o and its growth rate is 
limited by t he Ly apunov condition or Lindeberg condition ( Embrechts et aL[|l997t 
iBillingslevt Il995h . If the random variables are not independent, one can always 
group sets of such variables into independent and identical distributed clusters of 
variables and, again, the result would be a Gauss ian d i stribu tion. One exception 
is out from these Gaussian distributed situations (| Voitl l2003h : one where the size 
of these clusters is of the order of the total size of the system, i.e. when the second 
moment diverges, E(X 2 ) — > oo, violating Lindberg condition and, consequently, 
all the stronger ones. For such kind of dependent processes instead of a Gaussian 
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distribution, a a-stable Levy distribution with a < 2 and the typical heavy tails 
and power-laws is observed. 

On a more ontological approach, when modeling financial and economical 
networks, random variables are translated into agents. Agent-based models for 
describing and a ddressing the evolution o f markets has become an issue of in- 
creasing interest dFarmer and Folevl 20091) and appeals for further develo pments 
(jHaldane and Mavl l201lt I Johnson and LuxL l201l[ Samanidou et all 120071). The y 
enable one to access three important questions ( Bouchaud and Potters . 120031) . 
First, the system is able in this way to be decomposed into sellers and buyers, 
a common feature of all finance systems. Second, one does not need to assume 
the system to be in equilibrium, which indeed is what we observed in real stock 
markets. Third, by properly incorporating the ingredients of financial agents and 
the trades among them one can directly investigate the impact of trades in the 
price, according to some prescribe scheme. 

Here, we aim to choose a good model for the individual behavior of single 
financial agents, at a microscopic scale, in a way that the collective behavior gen- 
erates an output in accordance with the observed curves of macroscopic variables, 
namely the financial indices. 

Several of suc h bottom-up appro aches were thoroughly investigated. The Solo- 
mon-Levy model ( Lew et all 1 19941) defines each agent as a wealth function LUi(t) 
that cannot go bellow a floor level, given by U)i(t) > u>ou)(t) where w(i) is the 
agent average oj at instant t and wo is a proper constant. The imposition of the 
floor based on the mean field U)(t) means that on average (\wi(t) — u)(t)\) ~ N 
and, by basic statistics, var(cj(t)) ~ iV 2 . Consequently, the result of the Solomon- 
Levy model, despite the int eresting idea o f the introduction of a floor similar to 
what was done by Merton (jMertoij, Il973l) in the agent dynamics, will surely be 
a a-stable distribution with a power law heavy-tail, i.e., a < 2. A fo rtiori, perco- 
lation base d models like Cont-Bouchaud ( Cont and Bouchaudl . [20001) or Solomon- 
Weisbuch ( Solomon and Weisbuchl . Il999l) . by the nature of the phenomena, also 
brings up variations of system size order of magnitude and, obviously, will lead 
also to Levy-type distributions. If a multi-agent model does not give additional 
information and only provides the distribution a system size second moment, there 
is no particular surplus in the application of such a model compared with a direct 
fitting of an a-stable Levy distributi on. Moreover, all complex ity used to build in 
the individual behavior of the agents (|Lux and Marchesill2000t) will not contribute 
to the actual solution of the problem. 

Following the above considerations, we can eliminate from our reasoning every 
microscopic scenario that result on either weak or strong forms of CLT that will 
result on a Gaussian distribution. What remains? Every microscopic scenario that 
generates a system size order variance, with a a < 2 Levy stable distribution 
having the observed heavy tails. For everything else we can choose arbitrarily 
our model composition, which opens a rather broad panoply of possibilities. This 
mathematical freedom to choose the model of the underlying mechanisms is only 
constraint by the requirement of bridging the biological foundations of an economic 
system and the behavior of a market curve. With such bridging we are able to 
anchor information from the market curve to est ablished theoretica l princ iples of 
economy without additional assumptions like in ( da Cruz and Lindl . l201l|) . 

The question we are interested in is what is the minimum model for an economic 
system? In other words, what are the fundamental assumptions common to all 
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economic systems? Do they suffice for the emergence of statistical features of 
macroscopic variables such as heavy-tails? Next we argue that there are three 
fundamental assumptions for an economic system of agents. 

First, agents promote trading. Human beings are more efficient doing special- 
ized labor than being self-sufficient and for that they need to exchange labor. The 
more labor exchanges they make the more specialized they can be. The usage of 
the expression 'labor' can be regarded as excessive by economists, but readers can 
look at it as the fundamental quantity that is common to labor, money, wage, etc. 
that justifies the exchange. Something must be common to all these quantities; 
if not, we wouldn't exchange them. The physicists can regard such fundamental 
quantity as an 'economic energy'. 

Second, we only consume and produce a finite amount of the overall product. 
Thus, if an agent transposes that finite amount it loses his (consumption) trade 
connectivity in the economic environment in order to return to an admissible state 
that guarantee the finitude of the overall product within the system. 

Third, human agents are different and attract differently other agents to trade. 
As we see bellow, this difference should reflect some imitation, where agents tend 
to prefer to consume (resp. produce) from (resp. for) the agents with the largest 
number of consumers (resp. producers). The number of producer and consumer 
neighbors reflects supply and demand of its labor, respectively, and combining 
both kind of neighbors should suffice to quantify the price of the labor exchanged. 

In this paper we show that heavy-tailed return distributions are due to the 
economic organization emerging in a complex economic network of trades among 
agents governed under the above three assumptions. Further, we show that the 
power-law tails are characterized by an exponent that can be measured and is 
constrained by upper and lower bounds, which can be analytically deduced. The 
knowledge of such boundaries are of great importance for risk estimates: by de- 
riving upper and lower bounds, one avoids underestimates which enable the oc- 
currence of crisis unexpectedly and avoids overestimates which prevents profit 
maximization of the trading agents. 

We start in Section [2] by discussing how the biological foundations of the ex- 
istence of an economy are materialized into an economic organization topology. 
In section [3] we will show that an economic system behaves critically, transiting 
intermittently between two phases. Critical system combined with the underlying 
topology enable us to bound the behavior of the system, under topological con- 
straints. Such boundary behavior leads naturally to a maximum and minimum 
value on risk evaluation, improving the knowledge about the uncertainty of the 
market future evolution. Finally, we will show how empirical data and computa- 
tional results are in accordance with the analytical conclusions. 

2 System Topology in Attractive Fields 

If we look to corporate advertising or brand publicity actions one of the main 
characteristics is the usage of the volume of sales or the number of costumers 
as a motivating factor, though a large number of customers represent a higher 
demand and, consequently, a higher price, according to the Law of Demand and 
Supply. We feel encourage to buy something if we know beforehand that many 
other bought it already. The relative attractiveness of the demand is stronger 
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than the repulsive effect it has on the price. Som e auth ors in Economics call 



this phenomenon imitation ( Bouchaud and Potters! [2003), w hile physicist s call 



it preferential attachm ent and it was introduced by Simon ( Simonl Il955l) and 
developed by Barabasi ( Barabasi and Albertl Il999h . 



In this Section we argue that imitation or preferential attachment is an uni- 
versal principle in a multipoint attractive field environment, where each person, 
company or economic agent imposes an attractive field in the economical environ- 
ment and it can be expressed by a well known probabilistic law. 

For a better understanding, we take the human system as a general system 
of many bodies interacting through an attractive field. One starts by defining 
the potential function, based on a two body problem to get an Ham iltonian, the 
function that describes the dynamics of the system ( Greinerl 120031) . With such 



function and based on the existing interactions and their dynamic parameters - 
e.g. distance - the definition of the dynamics of the system is closed by means 
of principles like energy conservation and minimal action. In classical mechanics, 
what describes how planets orbit around stars or other planets? There is a well 
in the potential function that makes that one particular planet to be bonded to a 
particular start or bigger planet. But then two questions raise. 

First, why is a particular planet trapped in a particular potential well since in 
a many bodies problem there are many potential wells where the planet can be 
in? Second, what forms the potential well in a problem where everybody in the 
system attracts everybody else? 

If it is assumed that, besides the potential, no other kind of interaction or 
constraint exits, then the potential wells are formed exclusively by the same kind 
of interaction. Once the body is in a particular well, all related dynamics is well 
known since it is reduced to a specific two body problem. 

Suppose now that the two body problem is not known or it is not possible to 
formulate a general solution for a two body problem, like in the case of economic 
interactions. Additionally, suppose there exists a generic system of N — > oo bodies, 
not necessary in equilibrium, such that each body i exerts an attraction force Gi 
over the other ones, depending only on an intrinsic quantity Mi of itself and seen 
by the others. Thus, d = G{Mi). At each point A, the total set of bodies will 
produce a field which sums up the contributions from all agents: 

N 

g(A) = J2g(m>). (i) 

i = l 

When it is not possible to say more about G one expects that a body at A may 
orbit body i with a probability pi given by 

Expression ([2]) lays in the heart of our approach, since it reduces the many- 
body problem to a probabilistic expression, mainly because the solution for the two 
body problem is not known besides the assumption that there exists an attraction 
force between them. A first order approximation, G(Mi) ~ Mi, leads to 



Mj Mj 



±v±i ivi l 
P' ~ TvTT ~ — —■ ( 3 ) 
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Equation ([3j) incorporates some of the physics behind the two-body problem. Since 
Pi is dimensionless, the property Mi behind the local strength of the field is arbi- 
trary. 

Even not knowing which property Mi is, it is possible to measure the proba- 
bility pi from Eq. ([3j) . If Mi of a generic body i is relatively large when compared 
with the M's of the other bodies, then the number of bodies attracted by i will 
be large as well when compared with the other bodies. Assuming that the approx- 
imation in equation ([3]) is valid, then the relative number of bodies attracted by 
the body i represents the relative magnitude of Mi . As the number of bodies goes 
to infinity, the relative number of bodies effectively attracted by i converges to the 
probability for body i to attract other bodies. 

According to the Law of Large Numbers, this number of neighbors divided 
by the total amount of agents will converge to a theoretical probability that, 
in our particular case, is given by (|3]). Therefore with no reason for believing 
otherwise, the number of neighbors is taken as a relative measure of the field 
strength, Mi ~ k%, reducing Eq. ([3]) to 

"'TO 

The expression has been used a s assumption on several previous works with 
the name of "preferen t ial att achment" ( Barabasi and Albert . 19991 iDorogovtsev et all 
120001: iGonzalez et all 12006ft or "imitation" (jBouchaud and PottersL l2003h for the 
representation of the natural tendency that humans have to connect to the ones 
that already have more connections. Here we argued that it is a general rule for 
attractive interactions, at least on a first order approximation. 

There is an important consequence of this conclusion i n term s of the topology 
of systems, namely it is possible to state ( Barabasi et allll999l) that, if only this 
attractiveness is considered in the network, the probability of randomly finding a 
body with ki connections is 

P{k l ) oc kp (5) 

where 7 is o ne exponent, characte rizing of the system. We will show this result 
latter, as in (jBarabasi et all Il999h . since we will need it for the understanding of 
the criticality of the economic system. 

A power-law as the one in Eq. J5| is a rather general law for the topology 
(structure) of a system with different bodies out of e quilibrium, and it can be 
observed in many economic like syst ems as airports ( Li and Cail |2004|) or the 
World Wide Web (jBroder et all l2000h . 

Thus, the important conclusion from these section i s that, using the assumption 
that there is a biological surplus in exchanging labor ( Lipsev and Chrvstall lioOTj) 
- or representative forms like money - that generates the natural propensity of 
economic agents (people) to interact and that this agents are different from each 
other, the consequence is that the bodies will organize themselves in a topology 
given by Eq. ([5|). 



3 Non-Linearity and Critical Behavior 



It has been observed by some authors ([Sornetti [200I I Lux and Marches! lioooh 
the existence of critical behavior in stock market return distributions. Using a 
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non-linear threshol d of lev e rage s imilar to Merton's approach for the valuation 
of corporate debt ( Mertonl . Il973|) and the assumption that agents promote the 




exchange labor, here we show that the critical behavior exists in fact and it explains 
the heavy-tailed distributions of the stock market index curves. 

We define an economic connec- 
tion as an exchange of labor between 

two agents, say i and j, which dissi- „ pu- 

pates an amount of energy Ui . Agent J 9 ' 

% delivers an amount of labor Wij 
to agent j and gets a proportional 
amount of labor E%j = a%j Wij where 
aij is the 'exchange rate' of labor. 
Figure Q] illustrates the economic 

connection between two agents. Fig 1 iii ustration of CCO nomic connections. 

Being directed from i to j, agent i Details in text, 
takes the connection as a production 
or outgoing connection and agent j 

as a consumption or incoming connection. Assuming that each agent can connect to 
several other, since every agent will have the propensity to establish new economic 
links to get specialized labor from other agents, the "economic energy" balance 
equation on each agent i will be 

k iout k iin 
u *= £ W ij (l-a ij ) + J2 W m i(a m i-1) (6) 

where Lout is the set of neighboring agents to which i delivers labor and Lin the set 
of neighbors from which i gets labor. The energy in Eq. (JSj) has a straightforward 
parallel to the physical energy dissipated by the agent. Since energy and labor in 
this context have arbitrary units we approximate Wij = (W)t=o, representing Ui 
in units of (W) as 

Ui — ki^ ou t ki^i n -\- ^ ^ OLmi ^ ^ ^ij' (7) 

mEL in j€L out 

Assuming that the difference between any individual exchange rate of labor and 
the system average a = (otij) tot is relatively small, we can make a mean-field 
approximation as 

Ui = fi{ki,out &i,27i) (8) 

with ft = (1 — a). We must emphasize that we did not assume any additional 
assumption besides the ones that we need to justify the existence of an economy, 
we incorporate the economic or labor value in an energy function, since agents are 
physical entities that need to consume energy to generate labor. 

Since the energy consumption is necessary to establish economic connections, 
agents cannot leverage themselves to infinity, i.e. one agent cannot consume more 
then a certain amount of labor from other agents. Thus, we define a relative 
threshold of leverage U t h,i of an economic agent above which the other agents 
break their economic connections. We can regard this non-linear threshold as a 



s 



limit for default, as in credit risk modeling (|VasicekL[l987l) . Let us assume that 



this threshold is defined as a function of the total turnover defined as 

Ti — hi^out ~\~ (9) 

which measures the amount of trades, both consuming and producing, that in total 
an agent makes. The threshold must be defined relative to this economic impor- 
tance in a similar way as the u sage of financial ratios as a measure of probability 
of default ( Thomas et al l l2002l) . Thus each agent has a "relative" energy di which 



must obey 

di = I > d thti = ^. (10) 

The exchange rate of labor a.ij can be taken as a dimensionless price. Since 
from Eq. ([8]) one can look at ki n and k ou t as measures of supply and demand 
respectively, aij must grow monotonically with kout and drop monotonically with 
kin. Additionally, aij > and the first derivatives must decrease when \k ou t — 
kin\ — > 00. In other words, aij is a step function of (ki tOU t ~ fci,in) with aij to its 
maximal (resp. minimal) value for k ou t 3> ki n (resp. k ou t "C ki n ). 

So, our Ansatz for numeric simulation purposes is 

(11) 



13 1 _|_ e -(*W,i-fcta,j) 

which is basically a step funct ion with average value at Q max /2, similar to the one 
proposed by Cont-Bouchaud ( Cont and Bouchaudl . I2OO0I) . Henceforth we consider 



2 without loss of generality. 
When Eq. (|10[) is not fulfilled, the agent loses its incoming connections and will 
not be able to pay to its neighbors that will stop "working" for him. When this 
occurs the following is imposed to the collapsed agent i and its neighbors, indexed 
as j: 

U l ^rU l +fik liln (12a) 
Ti Ti - ki, in (12b) 
Uj -+Uj-P (12c) 
Tj Tj - 1. (12d) 



This "multi-agent" model of the economy is such that agents are not inde- 
pendent from each other. Their inter-dependence is here put in an intuitive form, 
through Eqs. (|12[) and basic assumptions in Economy justify completely Eq. (|10[) : 
(i) agents tend to deal with each other, (ii) it is not possible for the system to 
consume infinite energy and, for that reason, a limit leverage for each agent must 
exist and (iii) agents have an attractiveness forming an attractive field with a given 
topology, in which we focus next. 

Following the method by Barabasi- Albert- Jeong ( Barabasi et all Il999l) . at 
each time one new node enters the sy stem and new connec tions are formed by 
the imitation attachment mechanism ( Borland et all 120051) defined in Eq. (£[]). 
The growing rate of connections to each agent is given by 



(13) 
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with N the total number of agents in the system. Since each connection affects 
two agents, the denominator in Eq. (|13[) obeys k 3 = 2wit where £ is the total 

system age. Then = || with the solution 

k i {t)=m(j\ (14) 

where t% is the first instant when th e agen t i got into the economic network. Using 
the approach from ( Barabasi et all Il999f) . the probability for an agent to have k 
neighbors is 

d „„ , s d „, m 2 t, d (_ m 2 t \ 2m 2 t , _ a 



P(k) = < k) = -P(u >_.) = _ . k2{t + mo)J - mo + t ' 

(15) 

The exponent value three has an important meaning that we will address in 
Sec. [5] For now, we assume that the probability in (|15p has a general exponent 7, 
yielding P(k) ~ A; -7 . According to Eq. (|10p an agent will collapse, i.e. will lose 
their connections, if the following conditions are fulfilled: 



or more simply 



ki^out ki^in > Uth\ki^out ~\~ ki^in^) (l^V 
ki^out 1 ki^i n < Uth{ki,out 1 ~t~ &i,in) \^) 



< Ujki,in + 1 (18) 



with oj = ■ Meaning that the probability for node i to collapse is 



Pbr — P^L^ki^in *C ki^out < Uki,in + 1). (19) 



and ki t out are integers there is only one integer value for kt t out in the 
interval ki,i n < ki }0 ut < uki : i n + 1, given approximately by fci,out ~ wfei.in + \ 
and leading to 

P br « m(wfci, in ))" 7 , (20) 

since u)ki,i n ^> 1- 

Having a node that collapses, if its number of neighboring agents expected to 
collapse is bigger than one, the system would be in permanent destruction. On 
the other hand, if that expected number is lower than one, the system would grow 
indefinitely, since simultaneously agents tend to create more connections. In the 
first case one is lead to a scenario where agents are self-sufficient and in the second 
case one leaded to a non-realistic situation of infinite energy consumption. The only 
possible state of economy is the one where the expected value of neighbor agents 
to collapse is precisely one, which corresponds to a critical state. Obviously, this 
is a theoretical state. In practice, the economic system will travel from one phase 
to the other, growing and dropping as economists always said and as everybody 
always experiences for good or worse. 

Since the expected value must be equal to one then 

oo 

Yl k Hr> P ( k hn)Pbr = l (21) 
hi. =1 
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and from Eq. (fT5j) 

oo 

uP =m? ( fe ^n))" 27+1 ™ 7 C(27 - 1) (22) 

=1 

where £ is Riemann Zeta function. 

A similar derivation can be done for positive variations, instead of drops: in- 
stead of limiting the leverage to collapse, one limits the savings for a consumption 
chain of events that can lead to a heavy-tailed distribution with different exponents 
on the positive and negative side since the mechanisms are not dependent. 

With expression (122[) we show that, only by imposing first principles in Econ- 
omy, PI, P2 and P3 above, the economic system is trapped in a critical state 
defined through the relation between economic growth m, leverage level uj and the 
economic organization characterized by 7. Furthermore, the critical behavior leads 
to aggregate measures based on the exchange of labor, like returns, which have 
an a-stable Levy distribution. Differently from previous studies, here we focus on 
universal features of trading networks, without the need for more assumptions on 
agents like imposing them to be rational or irrational among other features. Nev- 
ertheless, we need to show that empirical distributions are in a ccordance with our 
concl usions. And, as we will describe next, they are in fact ( da Cruz and LindL 
l2010h . 



4 Exponent Amplification 



Since we know that the economical network is in a critical state, transiting between 
two phases we can look how the microscopic mechanisms that make neighbor 
agents to collapse generates the heavy-tails. Such mechanisms should be the same 
for the emergence of bubbles, but here we focus only in the negative part of the 
fluctuations distribution. 

The collapse of a node leads to 
the breaking of its consumption links 
or, equivalently, its neighbors' pro- 
duction links. Meaning that each 
collapse will provoke a chain re- 
action of size r, i.e., it originates 
a branching process as illustrated 
in Fig. [2] with a probabili ty ac- 
cordi ng to Otte r's theorem (lOtter , 
19491) . given bv (IOtterUl949t lHarrig . 




1996) 



1963:lAthreva and NevLfeoOl; Lau ritsen et al 



P(Z = r) cx r 



(23) 



Fig. 2 Illustration of a collapse chain reaction. 

where Z is the total number of nodes 
involved in a single branching process. Equation (123|) holds independently of the 
topology as long as the branching process is critical, but the number of collapsed 
agents in a real economic network is difficult to recount for. In the economical 
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context, one chain reaction is called a "crisis" and "avalanche" in the physical 
context. Henceforth, we address to them as crisis or avalanches indistinctly. 

What is measured when an avalanche occurs in such a real network is the 
number of links destroyed during the avalanche. This number of links accounts for 
a macroscopic property of the system, namely the overall product Ut which sums 
up all outgoing product of all agents. Therefore, we want to express P(Z) in terms 
of the total number of destroyed links. 
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Fig. 3 (a) Heavy-tails of stock indexes cumulative distributions and (b) the corresponding 
exponents. Distributions in (a) are shifted in the vertical axis for easy comparison. Solid line 
in (b) indicates the theoretical value obtained for 7 = 3 (see Eq. 1251 1. while dashed lines 
indicate the bounding values of m for 2 < 7 < 3. 



From the total number r of nodes included in one avalanche, since P(k) ~ k 1 , 
the number of nodes with kj connections involved in the avalanche is rij = rkj 1 

which corresponds to kjUj = rk~ 1+1 destroyed links in nodes with kj links. Thus, 
the total number of links destroyed is given by 



Kn 



-7+1 



(24) 
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where k m in and k m ax are the lowest and highest degree involved. Since the sum 
on the right side of (|24|) are the degrees, not the links, we can substitute each kj 
by ajKx, where aj are leverage dependent coefficient, which leads to r oc K^,, 

i.e. P{K T ) oc K~ = 7 . 

So, the fraction of avalanches of size — number of lost links - larger than size 
s is given by 

P(K T >s)oc x'^dx oc s -2 7+1 = s~ m (25) 



Equation (|25|) should hold for all types of trades and macroscopic observables 
of the economical product in one economic network. Figure [3] presents data from 
several stock market indexes. Figure [3^, shows the probability density functions 
(PDFs) of the logarithmic returns of each index, symbolized as x, where one can 
observe the heavy-t ails. While the hea vy-tailed shape of the return distributions 
was already known ( Mandelbrot!, Il997|) . at least within a good approximation, the 



explanation for their emergence was up to our knowledge not addressed till now. 

The result in Eq. f|25|) relates the microscopic economic relations - trades 
among agents - with the form of the heavy-tailed return distributions. This re- 
lation emerges from the interplay between the distribution of the returns, one 
macroscopic quantity chosen as the overall product Ut, and the topology of the 
social network. The economic connection was defined as one exchange of labor, in 
one of its forms. No additional assumptions were taken. Since we never made any 
conjectures regarding the specific type of economic interaction. 

Figure [3b shows the exponent m = ^7 — 1 in Eq. (|25|) for each one of the 
stock market indexes. All of them take values around the theoretical prediction 
m = 5/2 for a topology with 7 = 3. Further, all values lay between to bounds, 
rrimin = 2 < m < X = m m ax , as we explain in the next section. 



5 Distribution Boundaries 

In previous sections we saw how from first principles in Economy, namely that 
economic agents are impelled to exchange e conomic value in the form of labor, 



wealth, or other ( Lipsev and ChrvstalL 120071) and that they do that according to 



a imitation procedure, which we translated as a preferential attachment among 
agents, which yields the probabilistic law in Eq. Q. This was possible because it 
is not possible to know the nature of the labor Wi, and therefore Eq. Q contains all 
information about the labor of each agent necessary to establish the organization 
of the bodies in a general economic problem, where a financial problem is just a 
particular case. 

The form of the agent topology is defined through Eq. ([5|). But there is still an 
open issue: which values of 7 are admissible in an evolving network of economic 
trades and collapses? Observed economic netw orks show that typically 2 < 7 < 3 



( Albert and BarabasiL l20o3 iLi and CaiL l2004h Dut that empirical evidence is not 



enough to estimate the risk associated the exponent of crisis sizes. In this section 
we show the t opological causes that underlay the emergence of thos e bound values. 
We know ( Albert and BarabasiL [2003 iDorogovtsev et al.Ll2000h that economic 



like complex networks exhibit a self-similar topology that can be translated into the 
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probability power-law given in Eq. (J5j) • Because of this self-similarity, the same net- 
work topology is observed independently of the scale. The study of self-similarity 
in complex networks is done by getting the law that keeps nodes, links and their 
relationship invariant in scale. See illustrations in Fig. |U To this end, we separate 
the two cases of directed and undir ected networks an d use the renormalization 
method by Song, Havlin and Makse ( Song et all 12005^ . 

In directed networks, for each node links are either outgoing or incoming, the 
probability of a node to have k outgoing links — or correspondingly k incoming 
links - depends also in the scale p one is considering the entire network. The plot 
a) in Fig. [1] illustrates how unitary sets of agents compose a self-similar network 
of directed edges. 

Since in directed networks the number of probable states grows with the num- 
ber of nodes (see Fig. [3]) then we can derive the scale invariant law as 



dNP d { P , k) 
dp 







(26) 



meaning that the the node density in fc-sp ace NPd(p, k) is a scale invariant mea- 
sure, corroborating some previous results (|Song et all 120051 : l2006h . 

Plot b) in Fig. [3] illustrates the same as plot a) but this time for undirected 
edges. For undirected networks, however, the degree law Pb(p, k) is the probability 
for a node to have k start links or k end links indistinctly. Undirected networks 
can be regarded as compositions of two directed networks, each one with an equal 
directed degree law Pd(p,k). When one link is add or removed, affecting a node 
with degree k in one of the directed networks, it changes also a node with degree 
k'. 

As described in Fig. [3] we can see that the renormalization generates TV 2 ad- 
missible states when going from one scale to the next, since now each node can 



a) 




b) 




... 



Fig. 4 Illustration of renormalization in complex networks. Dashed connections are the ones 
that disappear as one goes down in scale, a) For a directed network, the number of probable 
states grows with the number N of renormalized nodes, since only the left nodes on the first 
step and the lower node on the second step contribute to the number of probable states, b) 
For an undirected network, the number of probable states grows with the square of the number 
of renormalized nodes N 2 , since all nodes in all steps contribute to the number of probable 
states. 



14 



connect to any other and simultaneously be connected from any other and the 
scale invariant law is 

dN2 ^ = 0. (27) 
dr 

In short, the self-similar transformation of the node degree, i.e. the number of 
links in a node will be ruled by 

N q P d {k)dk = N^P d (k p )dk p (28) 

where k and k p are, respectively, an original and renormalized quantity of links 
and AT and N p are the correspondent quantity of nodes. The exponent is q = 1 for 
directed networks and q = 2 for undirected networks. 

Since the power-law in Eq. §5$ is invariant under renormalization, 

P(k) -> P(k p ) oc k~ 7 (29) 



and if the renormalization of the links are expressed as ( Song et all 120051) k — > 
k p = a p k, renormalization will result in 

N" = N$dp +1 . (30) 

We now define Ib as a distance in the economic space between nodes, which 
equals the average number of links separating a randomly cho sen pair of agents . 
The fractal dimension ds of the network can the n be calculated ( Song et al.[ 2005) 
using a box-counting technique ( Falconet . 1 19971 ) as 



N p = Nig (31) 

and the links will scale as 

kp = kl g (32) 

which yields s p = l~j 3 dk . Thus, using Eqs. (J5j> , (|28|) . (|31|) . and (|32|) . the 7 value can 
be obtained as a function of the fractal dimension of the system as 

7 = 1 + 9^ (33) 



wher e q = 1 for directed networks and q = 2 for undirected networks ( Song et al.l 
l2005h . 

Next, we show that this result retrieves a topological constraint for the value of 
7 and, consequently, for the 'weight' of the heavy-tail in the crises size distributions 
given by Eq. ([5]). 

Let us suppose (wrongly ) that economic networks are unweighted, i.e. al l links 
are equal. It can be p roved (jAlbert and BarabasiL l2002t iDorogovtsev et all l2000l : 
iBarabasi et all Il999h that adding a node with a link to an existing node ch osen 
through a preferential attachment probability ( Barabasi and Albertl . Il999h will 



lead to a topology with degree distribution given by Eq. (J5j) with 7 = 3. That is 
coherent with our result (|33[) since adding one link for each node will make links 
to scale like nodes, meaning that ds = dfc and 7 = 3 because at the end of each 
link th ere is a node and for each node there is a link, except for very few seed 
nodes ( Albert and BarabasiL I2OO2I) . 
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If one assumes that connections between pair nodes do not depend on their 
connectivity or, equivalently, all nodes have the same connectivity then for a set 
of N nodes we find N 2 — N links since all nodes connect between them. Thus, the 
links scale with N 2 — N ~ N 2 when N — > oo, yielding dk = 2ds, i.e. 7 = 2. 

Since in an economic network links are not all equal, the considerations above 
imply that we should expect values of 7 < 3 in real economies. Preferential attach- 
ment and maximum number of random links are the two extreme situations of a 
spectrum of possibilities corresponding to values of 7 in the interval 2 < 7 < 3. 

Additionally, if we assume that more than one type of connections can coex- 
ist, the final 7 can assume a value under the limit 3, which supports studies on 
empirical economic indices ( da Cruz and LindL l2010h . 

This is an important result since, independently from the complexity associated 
with economic processes and with the network of economic relations, we can bound 
the return distribution and, thus, find a closed risk measure. 



6 Conclusion 

In this paper we establish the bridge between basic economic principles and the 
heavy-tails of the return distributions whose exponent values are bounded. The 
boundary values measure the model risk associated with any model. 

We showed that using only the assumption that economic agents attract each 
other to exchan ge labor, which is the fund amental assumption for the existence 
of an economy ( Lipsev and ChrvstalL l2007f ). then economic links tend to be dis- 



tributed according to a power-law. 

We argued that return distributions with heavy-tails observed empirically, 
emerge from the aggregation of dependent random variables, which has variance 
diverging with the system size. Consequently, for preventing the system from di- 
verging, there is a wealth threshold for e ach economic agent, exactly as predicted 
by Merton for corporate debt modeling ( Mertonl Il973l ). 



Combining both findings - wealth threshold and power-law - we showed that 
the exponent is limited to 2 < 7 < 3. Since the exponent m of the heavy-tails for 
cumulative return distributions (sizes of crises) is given by 37/2 — 1 , we showed that 
2 < m < 3.5. Consequently, the heavy-tails observed in the return distributions 
assume values between 3 and 4.5. 

These findings help to solve the controversy about Mandelbrodt hypothesis 
(|MandelbrotL Il997h that the distribution of financial returns are explained by 



Levy distributions, and therefore would yield an exp onent smaller than 3 . Some 
arguments against Mandelbrot hypothesis were based ( Borland et al.l [20051 in em- 



pirical measures of the return distributions, yielding exponents larger than 3. As 
summarized in the previous paragraph, our findings show the same range of expo- 
nent values. Still, we think that Mandelbrot was right. Though what we measure 
for the time series of returns are links between agents (corresponding to exponents 
larger than 3), the random variables behind Levy distributions are the agents, 
which corresponding to the exponent 7 (smaller than 3). This work bridges these 
two apparent contradictory positions. 

Finally, since the exponent is bounded, the total risk associated with the pro- 
cess being observed is also bounded between a lower and one upper boundary 
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values. This is perhaps the most interesting result from our work: Having bound- 
aries for the heavy-tailed return distributions we can actually measure model risk, 
since the economic organization has topological limits. 
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